Abrupt appearance of the domain pattern and fatigue of thin ferroelectric films 
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We study the domain structure in ferroelectric thin films with a 'passive' layer (material with 
damaged ferroelectric properties) at the interface between the film and electrodes within a continuous 
medium approximation. An abrupt transition from a monodomain to a polydomain state has been 
found with the increase of the 'passive' layer thickness d. The domain width changes very quickly 
at the transition (exponentially with d~ 2 ). We have estimated the dielectric response dP/dE (the 
slope of the hysteresis loop) in the 'fatigued' multidomain state and found that it is in agreement 
with experiment, assuming realistic parameters of the layer. We derive a simple universal relation 
for the dielectric response, which scales as 1/d, involving only the properties of the passive layer. 
This relation qualitatively reproduces the evolution of the hysteresis loop in fatigued samples and 
it could be tested with controlled experiments. It is expected that the coercive field should increase 
with decreasing lateral size of the film. We believe that specific properties of the domain structure 
under bias voltage in ferroelectrics with a passive layer can resolve the long-standing 'paradox of 
the coercive field'. 
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Recent studies of thin ferroelectric films have revived 
interest in properties of the ferroelectric domain struc- 
tures. It became clear that some basic aspects of these 
properties remain unexplored, whereas they are of key 
importance for applications. For instance, the progres- 
sive loss of switchable polarization after repeated switch- 
ing cycles, i.e. polarization fatigue, is a serious problem 
in device applications of ferroelectrics. In spite of exten- 
sive studies, the physics of the fatigue remains poorly un- 
derstood. Various mechanisms were considered over the 
years |l]-||] • In many cases the deterioration of the switch- 
ing behavior, like the loss of the coercive force and of 
the squareness of hysteresis loop, were attributed to the 
growth of a 'passive layer' at the ferroelectric-electrode 
interface [l|^||, or to the pinning of domain walls by 
defects [0,4]. In this paper, we study the effects of a 
passive layer (material with damaged ferroelectric prop- 
erties) . We assume an ideal infinite ferroelectric film and 
treat it within the continuous medium approximation. 
The model exhibits very interesting properties relevant 
for real systems. 

The passive layer leads to the appearance of a depolar- 
izing field in the ferroelectric, so that the system would 
tend to transform to a polydomain state in order to re- 
duce the field. There is, therefore, a transition between 
the monodomain state, when the thickness d of the pas- 
sive layer is zero, to a polydomain state otherwise. What 
is the nature of this transition? This is the main theo- 
retical question we address in this paper. 

For a short-circuited infinite ferroelectric plate, as we 



shall show, the domain structure exists for any value of 
the passive layer thickness, in disagreement with some 
earlier results ||. A surprising feature of the present 
study is that the formation of the polydomain state sets 
in abruptly with the appearance of the passive layer. The 
period of the domain structure tends to infinity exponen- 
tially when the passive layer thickness d goes to zero. The 
theory allows one to study the evolution of the domain 
structure in the asymptotic regime of very wide domains 
and calculate its dielectric response. 

The slope of the resulting P — E (polarization - elec- 
tric field) hysteresis curve at E — is directly related 
to the thickness d and dielectric constant e g of the non- 
ferroelectric layer. We show that there is a universal re- 



lation (14) dP/dE cc e g /d. This relation is in good corre- 



spondence with available experimental data. It is impor- 
tant that even without pinning by defects, the response 
of this structure to an external bias voltage is rigid. The 
implication for real systems is that with the growth of the 
passive layer the hysteresis loop very quickly deteriorates 
and looses its squareness. 

The properties of the domain structure in the pres- 
ence of a passive layer are important for the problem of 
switching in ferroelectric thin films. If we were to apply 
the external field in the opposite direction to the spon- 
taneous polarization in a monodomain state, the equi- 
librium size of the domain with the 'wrong' polarization 
would eventually become less than the film size. Then a 
new domain wall would appear in the film and we get a 
switching. Obviously, the coercive field for this mecha- 
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nism should then increase with decreasing lateral size of 
the film. 

We shall consider a ferroelectric film under a bias volt- 
age U with thickness I separated from the top and bottom 
electrodes by passive layers with thickness d/2 (Fig. 1, in- 
set). For a system with a given voltage on the electrodes, 
the equilibrium distribution of charges and polarization 
corresponds to the extremum of the thermodynamic po- 
tential F, which includes explicitly the work produced 
by the voltage source(s) and the energy of the electric 

field easiai: 



F = F, 



LGD 



dVE 2 /(8Tr) 



electrodes 
a 



(i) 



Here -Flgd is the standard Landau-Ginzburg-Devonshirc 
(LGD) functional at zero electric field |Io| , ^2| , E is the 
electric field, while e a (<p a ) is the charge (potential) on 
the electrode a. The last term in Eq. ([!]) accounts for the 
work of the external voltage source(s). 

We are interested in a case when the ferroelectric has a 
spontaneous polarization P s \\ z (Fig. 1, inset). The film 
has a dielectric constant e c (e a ) in the z-direction (in the 
xy-plane), and the dielectric constant of the passive layer 
is e g . We select the x-axis perpendicular to the domain 
walls. The potential ip (E = ~Vp) satisfies the following 
equations in the ferroelectric and the passive layer, 



dx 2 



dz 2 



0. 



°<Pg 



dx 2 dz 2 



0, 



(2) 



with the boundary conditions p = — (+)U/2, z 
+(-)(« + d)/2, and 



iff = ip g , at z = 1/2, (3) 



where the subscript / (g) denotes the ferroelectric (pas- 
sive layer, or a vacuum gap). Here a is the density of 
the bound charge due to spontaneous polarization at the 
ferroelectric-passive layer interface, a = P ns = ±P S , de- 
pending on the normal direction of the polarization at the 
interface, alternating from domain to domain, Fig. 1 (in- 
set). We have assumed a usual separation of linear and 
spontaneous polarization, so that the displacement vec- 
tor is Di — €ikEk + P S i, where i,k = x,y,z, and the 
dielectric response in uniaxial, Fig. 1 (inset). 

We shall assume that the film is split into domains 
with alternating z-component of polarization and widths 
a\ and ai so that the pattern is periodic with a period 
T = a% + a%. This assumption is perfectly reasonable, 
since the bound charges at the boundary between the 
ferroelectric and the passive layer are fully compensated 
by the image charges at the electrodes. In this situation 
the domains would not branch (see e.g. |12), p. 219). The 
solution of the Poisson equations (|^) is then readily found 
by Fourier transformation 
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FIG. 1. The domain width a in a ferroelectric capacitor 
versus the passive layer thickness d for different separations 
between the electrode plates L, e g — e a — 10, e c = 200 at zero 
bias. Dashed lines show the behavior of the Kittel width ajf. 
All quantities are scaled to W, the domain wall thickness 
[12]. Inset shows schematics of the electroded ferroelectric 
film (capacitor) under bias voltage U. Note that the transition 
to very wide domains at small d is very sharp, even when 
shown on a logarithmic scale. 
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where n is the number of the Fourier component, and 
a = /, g. With the assumption that P s = const inside 
the domains, the free energy Eq. |l]) becomes F = Fq + 
\ /fe dAaLp ~i]J2a e afa, where F = F LGD [P S ] is the 
free energy of the ferroelectric with the polarization P s 
which contains the surface energy of the domain walls, 
and the integral is taken over the interface between the 
ferroelectric and the passive layer. It is useful to treat 
the homogeneous part (n = 0) of the potential and the 
surface charge separately from other terms. Then the 
result for the (dimensionless) free energy per unit volume 
is 



/ 



F 



S (4n5d/L) - e g u 
2 



ALP 2 2 (e c d/0 
e g u (e c uL/l) + 4ir5 
8tt (e c d/l) + e g 



+ fK(t) + f x (t,S,u), (8) 



where the dimensionless bias voltage is u = U/P S L, and 
the first two 'homogeneous' terms describe the growth 
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of the domains with the 'right' polarization (along the 
external field) at the expense of domains with opposite 
polarization, is the Kittel-like free energy [ p^|JTl| 

defining the dimensionless period of the domain struc- 
ture t = (ai + Oa)/L, and f x (t, 5, u) is the cross-term Q 
describing, in addition to the first two terms, the appear- 
ance of a net polarization in the film under bias, while 
L = l+d is the distance between the electrodes, Fig. l(in- 
set). To find the equilibrium distribution of charges and 
polarization, this functional should be varied with respect 
to 5 and t at constant bias voltage (parameter it). 

The domain structure with the passive layer can be 
determined simply from a minimum of Jk in the case of 
zero bias voltage, when a\ = a-i = a, and 



Lf K = 



Al 32a 



E 



(9) 



e g coth 



(2j+lM 
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^coth 



(t) 



1/2 



(2j+l)^ 
2a 



For a free standing film, d > I, Eq. (^|) has a standard 
(generalized Kittel) solution 



I 1/2 



a = a K = [ir 2 (e g + AZ/28C(3)] ' 

= [0.3(e ff + v^)AZ] 1/2 , (10) 

where A is the (temperature dependent) characteristic 
length, which is proportional to the product of the inverse 
of the dielectric susceptibility \ in the FE phase and the 
domain wall thickness W pT| , p2| . 

Most interesting is the ferroelectric film with a narrow 

passive layer, d <C I. When ( |& J |^ > 1, the second 

coth in (^) can be replaced by unity and neglected alto- 
gether in all terms in the sum with j < a/(jrd), since the 
first factor in the denominator is e g coth [?p(2j + 1)] ~ 
% ^> 1. We see that the factor a cancels out from all 
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those terms in the sum. The dependence of the second 
term (the so-called energy of emergence per unit area) 
on a becomes very weak, and the free energy minimum 
would correspond to very large domain widths, a 3> a^ , 
Fig. 1. The domain size may exceed the lateral size 
of the sample, and then we would have a monodomain 
film. This surprising behavior is the result of the effec- 
tive screening of the bound charges on the ferroelectric- 
passive layer interface by the metallic electrodes |l5| . 

To analyze the transformation between monodomain 
and multidomain states, we first consider an exactly solv- 
able model with e„ = »/ 



in the regime l/a> 1. Then 



Lf K = Al/a 



32a 



E 



Al/a H 1 



U {2j + 1)3 coth 
where 



(2j + l)7Td 

2a 



,__ e A(2j+l)77d/a 2 t p 



up to terms ~ where ^ = ird/a -C 1, and p = (3 + 
ln4)/2 [Q. Minimizing the free energy fx, we finally 
obtain the width of the domains in the case of a narrow 
passive layer 



ircl 
2?/2 



exp 



(12) 



This unusual solution is extremely sensitive to the thick- 
ness d of the passive layer, so that the domain width 
increases explosively for narrowing passive layer. This 
conclusion is general, it applies at a > I too, since the 
cancelation of a persists in this case as well, as illustrated 
by the results of the exact calculation in Fig. 1. Note 
that a similar case of a ferroelectric capacitor with non- 
ferroelectric insulating layer at the electrode surface has 
been considered recently fT^ | , but the exponential growth 
of domain thickness ([Ej) was overlooked, and only the 
standard Kittel expression was obtained. 

We have seen that the ferroelectric film splits into do- 
mains as a result of the growth of a passive layer at the 
electrode surfaces. This substantially decreases the coer- 
cive field. It will be defined by only the strength of the 
domain wall pinning. We can assess the validity of this 
scenario of fatigued state by considering an ideal case 
with no pinning and comparing the estimated slope of 
the hysteresis loop with experiment . To this end, wc shall 
estimate the value of dP a /dE, where P a is the apparent 
(measured) polarization (surface charge) P a = o~q = P S S, 
Eq. (g), and E is the electric field (in fact, it is more 
convenient to find dP a /dU = L~ l dP a /dE). The linear 
response to bias voltage is obtained from Eq. (||) , which, 
for small bias U, takes the form 



f = fo(t) + -S6 2 - RuS, 



(13) 



where S is the domain structure stiffness with respect to 
changes in apparent polarization P a , and the coefficient 
R defines the response to the bias voltage. The slope of 
the hysteresis loop L(dP a /dU) — R/S is shown in Fig. 
2. Since the cross-term f x in Eq. (8) changes only the 
numbers, we can obtain a simple analytical approxima- 
tion by neglecting this term to the first approximation. 
The result is 



dPa 
dU 



And' 



(14) 



7T£ 



This remarkably simple relation is universal, and de- 
pends only on properties of the passive layer [the ex- 
pression reads eoe g /d in SI units]. It is also rather ac- 
curate in comparison with the exact numerical calcula- 
tion, Fig. 2. It is well known from experiment that, as 
fatigue progresses (d grows), the slope of the hysteresis 
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FIG. 2. The calculated slope of the hysteresis loop, 
W(dP a /dU), in a fatigued sample in comparison with the 
approximate expression.Eq. (0). The parameters of the film 
are the same as in Fig. hi Insets show schematics of the very 
fast decrease of the slope of the hysteresis P a — E loop with 
the growth of the passive layer (fatigue). 



loop steadily decreases H||l8|. This behavior follows 
immediately from our Eq. (14), the slope of hysteresis 
loop decreases very rapidly with increasing d, Fig. 2. 
From experiment fl8| , we can determine the slope in PZT 
films as dP a /dU ~ 15/iCcm V -1 . One can estimate the 
corresponding thickness of the passive layer from Fig. 
2. Equation ( |l4| ) gives for the passive layer thickness 
d ~ 2 — 7 nm with e g = 3 — 10, the thickness of the 
film was L = 256 nm p8[ . This d may be somewhat un- 
derestimated since (i) the experimental slope is affected 
by pinning of the domain wall, which we neglected, and 
(ii) Eq. ( |l~4] ) somewhat underestimates the exact result 
(Fig. 2). The result shows that already a tiny passive 
layer has a profound effect on the hysteresis loop and its 
slope rapidly diminishes with the growth of the passive 
layer. Thus, in thin films with small lateral dimensions 
the transition into a fatigued state with domains will be 
discontinuous and the growth of the passive layer can 
easily trigger this transformation. 

Within our model there is no hysteresis for an infinite 
film and the dielectric response describes the slope of the 
hysteresis loop. It is important to emphasize that the 
film becomes more 'dielectrically rigid' because of the 
growth of the passive layer alone, even without taking 
into account either bulk or surface pinning of the domain 
walls. The estimated slope of the fatigued system is in 
good correspondence with experimental data, so that the 
growth of the passive layer might be the main origin of 



fatigue. 

We believe also that the specific features of properties 
of the domain structure in the presence of the passive 
layer may be a key to resolving the famous 'paradox' of 
the coercive field. All the estimates for creating a nucleus 
with reversed polarization give extremely large values of 
the energy barrier, E* ~ 10 7 k B T |||-§l|. Such a huge 
barrier makes a thermal nucleation impossible, and this 
makes the related estimates of the coercive field question- 
able 

In our scenario, however, the coercive field is defined 
by the energy barrier for the domain wall to 'enter' the 
sample. It is likely to be much smaller than the classic 
barrier for nucleation E* , and this possibly resolves this 
long-standing paradox. We hope that these new ideas 
about the problem of switching will prove fruitful. For 
instance, we expect that in thin films the coercive field 
will be increasing with decreasing lateral size of the sam- 
ple, since the size of a domain with a 'wrong' polariza- 
tion should become smaller than that size. This should 
be contrasted with the usual conclusion that the coercive 
field grows with decreasing thickness of the sample [ p2[ , 
irrespective of its lateral size. It would be very interest- 
ing to test this prediction experimentally, especially in 
view of ever smaller sizes of films used in devices. 

The effect of the passive layer on switching in the ferro- 
electric films is a specific case of a more general problem: 
the effect of inhomogeneities in the system on its ferro- 
electric properties. As we have shown elsewhere 23 1, a 
similar mechanism produces the fatigue of a film without 
any actual passive layer, in presence of a depletion charge 



in ferroelectric-semiconducting films 24 1. 

It is worth emphasizing that none of our conclusions 
depends on a particular form for the domain structure 
(which may be of a checkerboard type, etc.) This is be- 
cause the effect of partial screening of the depolarizing 
field by electrodes obviously applies to any domain pat- 
tern. 
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hospitality of Quantum Structures Research Initiative at 
Hewlett-Packard Labs (Palo Alto) during the course of 
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